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It is well-known that in thin ferromagnetic film with a net magnetization perpendicular to the film 
the collinear arrangement of spins is unstable in an in-plane field H smaller than its saturation value 
H c . Existence of a stripe phase was proposed with elongated domains of alternating direction of 
magnetization component perpendicular to the film. We consider in the present paper the strong- 
field regime H < H c and discuss the minimal microscopic model describing this phenomenon, 
two-dimensional Heisenberg ferromagnet with strong easy-axis anisotropy and dipolar forces. The 
noncollinear (stripe) phase is discussed using the technique of the Bose-Einstein condensation of 
magnons. Some previously unknown results are observed concerning the stripe phase. Evolution 
is established of the spin arrangement in the noncollinear phase upon the field rotation within the 
plane. We find a rapid decreasing of the period of the stripe structure as the field decreases. We 
demonstrate that spins components perpendicular to the film form a sinusoid in the noncollinear 
phase at H m H c that transforms to a step-like profile upon the field decreasing so that the domain 
wall density decreases from unity to a value much smaller than unity. The spin-wave spectrum in 
the noncollinear phase is discussed. 

PACS numbers: 75.70.Ak, 75.30.Ds, 75.10.Jm, 75.10.Dg 



I. INTRODUCTION 



Numerous fascinating magnetic properties of ultrathin films and (quasi-) two-dimensional (2D) magnetic materials 
attract much attention now that is stimulated also by the technological significance of such materials ^ One of the 
most intriguing experimental findings was the discovery of the spin reorientation transition (SRT) in thin ferromagnetic 
films. 3 Recent experimental and theoretical investigations have shown that SRT can be induced by the temperature, the 
film thickness or by an applied magnetic field. In most cases the rotation occurs during SRT from the perpendicular 
to an in-plane direction (or vice versa) of the film magnetization. In many materials with the net magnetization 
perpendicular to the film this reorientation is accompanied by stripe phases with elongated domains of alternating 
magnetization direction which are separated by domain walls with a finite width. Such transition was obtained first 
in Refi^ in Co/Au(lll) films. The reader is refered to comprehensive reviews^ of numerous recent experimental and 
theoretical works in this field. It is well established now that the origin of SRT and stripe phases is in the subtle 
interplay between the short-range exchange, anisotropy and long-range dipolar interactions between spinsi^ Then, a 
realistic theoretical model of low-dimensional magnetic systems must include these three kinds of interactions so that 
the Hamiltonian has the form of the minimal microscopic model showing SRT and the stripe phase 
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where we direct y axis perpendicularly to the lattice as it is shown in Fig. [T] and the last term describes the Zeeman 
energy in the in-plane field discussed below. 

It is demonstrated in Ref. 5 that at H — and T = the stripe arrangement of spins in the model ([1]) has slightly 
smaller energy than any collinear arrangement if 2 A > qwq (i.e., if the easy direction is normal to the film), where 



3t>o 

a = 



E4 (3) 



8tt ^ Rf 

is a constant that is equal approximately to 1.078 for the simple square lattice, 

uj = ^(g^i) 2 (4) 



is the characteristic dipolar energy and we set the lattice spacing to be equal to unity. Meantime it was found in Ref£ 
that the period of the stripe structure rises rapidly as the value of the anisotropy increases so that the domain width 
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FIG. 1: (Color online). Two phases of the model |T} in strong magnetic field H: those with collinear and noncollinear spin 
orderings at H > H c and H < H c , respectively. Off- plane spins components are shown in red in the stripe phase (i.e., at 
H < H c ). It is obtained in the present paper that these components form a sinusoid at H ~ H c that transforms to a step- like 
profile upon the field decreasing as it is shown in Fig. [4] Period of the stripe pattern rises as the field decreases as it is presented 
in Fig. [S] 

is macroscopically large (i.e., larger than 10 8 lattice spacing) for 2A > lAaujQ. It is implied hereafter that J ^> ujq, 
as it usually is. We assume in the present paper that 2 A > lAatoo and consider macroscopically large unidomain 
sample with all spins directed perpendicularly to the film at H = 0. It is assumed below also that J A. 

Existence of a stripe phase in this case was pointed out long time age 6 in an interval of in-plane magnetic field 
H' < H < H r , where H r oc A is the saturation field, 



-Dk 2 



(5) 



D is the spin- wave stiffness that is equal, in particular, to SJ for simple square lattice with exchange coupling constant 
J between neighboring spins, and 



Sujq 
4D ' 



(6) 



The experimental evidence for the in-plane field induced collinear phase instability has been reported recently in 
Ni(001) ultrathin films at room temperature.— It was pointed out in Ref.— that the spin-wave spectrum calculated 
in the linear spin- wave approximation assuming a collinear spin arrangement is unstable when H' c < H < H c . 
Importantly, the instability takes place not at zero momentum but at a finite incommensurate one k c o which direction 
is not fixed in the linear spin- wave approximation and which value is given by Eq. ([6]). As a result it was proposed 
that the model ([T]) shows three phases in a nonzero in-plane field, collinear one at H > H c , collinear canted phase at 
H < H' c , and the noncollinear one appearing at H' c < H < H c . The properties of the noncollinear phase have not 
been clarified yet. It was proposed in Ref£ that there should be a domain pattern at H' c < H < H c with the period 
equal to 27r/fc c o forming by spins components perpendicular to the film similar to that discussed in RefJi for H = 0. 

It should be pointed out however that the stability of the spin-wave spectrum is not a criterion for the corresponding 
state to be the ground state. Then, it is not clear at which field H' c the transition takes place from the stripe phase 
to a canted collinear one. It is not clear also whether such a transition happens at all because the stripe phase has 
lower energy at H — 0, as it is shown in Ref£. 

Despite great recent interest to exotic non-collinear phases, that of the model JTJ) in strong-field regime has not 
been discussed thoroughly yet. The aim of the present paper is to fill up somewhat this gap. We discuss first 1/S 
corrections to the spectrum in the collinear phase at H > H c and show that they fix the direction of k c o at which 
the spectrum instability takes place and which absolute value is given by Eq. ((6|). We find the dependence of k c o 
direction on the field direction within the plane. 

The noncollinear phase is discussed below using the technique of the Bose-Einstein condensation of magnons sug- 
gested in Ref.-. The transition to the stripe phase at H = H c corresponds within this technique to a "condensation" 
of magnons at states characterized by a momentum k c ||k c o and its harmonics (i.e., nk c , where n is integer). It is 
the unusual result of the present analysis that the vector k c varies with the field (remaining parallel to k c o) so that 
k c = k c o at H = H c and k c < k c o at H < H c . This unusual behavior is related to condensation of magnons at states 
characterized by harmonics of k c . Remember that the condensation in antiferromagnetio 8 - and ferromagnetic 9 systems 
takes place at antiferromagnetic and zero (commensurate) vectors, respectively, which do not depend on magnetic 
field. Thus, we find that the period of the stripe structure given by 2n/k c rises very fast as the field decreases. 
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We demonstrate that the energy of the stripe phase is lower than the energy of any collinear spin arrangement at 
H < H c but the difference between these energies decreases as the field reduces. 

We find values of spins components perpendicular to the film forming the domain pattern, ft is shown that spins 
components perpendicular to the film form a sinusoid in the noncollinear phase at H ss H c that transforms to a step- 
like profile upon the field decreasing so that the domain wall density decreases from unity to a value much smaller 
than unity. 

We discuss below also the spin-wave spectrum in the noncollinear phase. 

The rest of the paper is organized as follows. We discuss the Hamiltonian transformation and technique in SecHTl 
Magnon spectrum is discussed in the linear spin- wave approximation in Sec. lIIIl First 1/5 corrections to the spectrum 
at H > H c is discussed in Sec. IIVI Properties of the noncollinear phase are considered in Sec. [V] Sec. I VII contains 
our conclusion. One appendix is added with some details of calculations. 



II. HAMILTONIAN TRANSFORMATION AND TECHNIQUE 



Qk = u o ( -5 P p - t^-^t ) i where Pifi = x > *> (8) 



It is implied below that the in-plane magnetic field is directed arbitrary relative to the lattice. Taking the Fourier 
transformation we have from Eq. {1} 

H = -\ E + Qu) S A - aJ2s v A - hVks* , (7) 

k k 

where J k = Ji m exp(ikR; m ), — Qff n exp(ikR; m ) and VI is the number of spins in the lattice. Dipolar 
tensor Q£r possesses the well-known properties'^ at k <C 1, which are independent of the lattice type and the 
orientation of x and z axes relative to the lattice, 

p0 fa k k p k(3 

Qjf = u ^-^a+^ Syp, where (3 = x,y,z, (9) 

where a and wo are given by Eqs. ([3]) and ([4|, respectively. It is seen from Eqs. (fTJ), JSj) and ([9]) that dipolar forces 
lead to easy-plane anisotropy in the energy of the classical 2D FM with y to be a hard axis^S Quantum and thermal 
fluctuations lead also to an in-plain anisotropy (order- by-disorder effect)] 11 ! 12 ' 13 

We use in the present paper the well-known representation of spins components via Bose-operators a and a^'J^ 

Sr*VMat(l-&Y S|^|l-|J a „ Sl = S-4*i (10) 

which can be obtained from the Holstein-Primakoff transformation^ by expanding the square roots up to the first 
terms. 

To perform the calculations it is convenient to introduce the following retarded Green's functions: G(o>,k) = 
(ak,a k )^, F(u,\l) = (a k ,a_ k ) w , G(a>,k) = (al k ,a_ k ) w = G*(-w,-k) and F^(oj,k) = (al k ,a£)« = F*(-oj,-k). 
We have the following set of Dyson equations: 

G(w,k) = G(°)(«,k) +G( )(a;,k)E(«,k)G(w,k) + G(°V,k)[B k + n(w, k)]Ft( w , k), 

Ft(w, k) = G {Q \oj, k)E(w, k)Ft(w, k) + G (0) (w, k)[S k + H^w, k)]G(w, k), (H) 

where G( )(a;, k) = (w — i? k + i5) _1 is the bare Green's function, E(w, k) = E(— oj, — k)*, II (a; , k) = IP (-cj, — k)* are 
the self-energy parts and i? k and £? k are coefficients in the bilinear part of the Hamiltonian 



Ti-2 = ^2 ^ka k a k + (a k a_ k + a k al k ) 
k <- 

which are presented below. Solving Eqs. (jTTJ) one obtains 

u + E k + E(w, k) 



(12) 



G(w,k) = 



X>(w,k) 



i? k + n^,k) 

F( "' k) = %kj ' (13) 
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where 

V(fi),k) = (uj + i5) 2 -el -Jl(w,k), (14) 
4 = K- Bl L ^ _ (15) 

n(u, k) = £ k (E + s) - s k (n + n*) - (w + - s) - nnt + ss, (16) 

and e k is the spin- wave spectrum in the linear spin- wave approximation. Quantity Q(lu, k) given by Eq. (|16[) describes 
renormalization of the spin-wave spectrum square. We find fl(uj, k) within the first order of 1/5 in Sec. IIVI 



III. CLASSICAL MAGNON SPECTRA 

A. H > H c 

Let us assume first that the field is so strong that all spins lie within the plane (see Fig. [T]). After substitution of 
Eqs. PH|) into Eq. ([7]) the Hamiltonian has the form H = E + Yli=i %i wnere E Q is the ground state energy and Hi 
denote terms containing products of i operators a and . Tii = because it contains only Q^f with p ^ (3. For Ti.2 
one has Eq. (112j) . where 

s k = 5(j -j k )-f (or+of - ? f^)+ ff -^( 1 -^) 

^ Dfc 2 - ^cos 2 ^ k + fl" - A5 (l - A) + (17) 

where k is the angle between k and the magnetization and the expressions after « are approximate values of the 
corresponding quantities at k <C 1. We find for the square of the magnon spectrum from Eqs. (fT5j) . (jTTJ) and (fTg|) in 
accordance with the well-known result^ 

^k, (19) 

(20) 
(21) 

where Hc^ is the classical value of the saturation field and k c0 is given by Eq. ©. The second term in the first 
brackets in Eq. (fT9")l is negligible in the considered case of A ^> cjo- The last term in Eq. (TIT))) is formally of the next 
order in 1/5 compared to the first one. We show below that corrections from the Hartree-Fock diagram presented 
in Fig. HJa) cancel it and leads to renormalization of the critical field He ■ We point out also that the anisotropy 
constant A should be accompanied by the factor 1 — 1/25 in all expressions for observable quantities because there 
is no one-ion anisotropy for 5 = 1/2M- It is seen that in the zeroth order in 1/5 (the last term in Eq. ([IT)]) should 

be discarded) the second bracket in Eq. (flT))) is positive at H > Hi and the spectrum is stable. On the other hand 

at H < Hc°^ the instability of the spectrum arises at momentum k c o. Notice that the direction of k c o is not fixed 
within the spin-wave approximation. Meantime we show below that there are anisotropic corrections to the spectrum 
of the first order in 1/5, which fix the direction of the vector k c0 at which the instability arises leaving, however, the 
equivalence between k c o and — k C Q. 



4 (V + ^fcsin 2 k + l/j (z?(fc-fc c0 ) 2 +A (0) )- 

A<°> = H-Hf\ 

H<® = 2SA(l-^-)-Su;oa + Dk 2 c0 , 
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b) 




FIG. 2: Diagrams of the first order in 1/S for self-energy parts. Diagram (a) stems from four-magnon terms in the Hamiltonian 
whereas (b) comes from three-magnon terms. 



B. 



H<H C 



Simple calculation leads to the following result for the spectrum at k -C 1 assuming the collinear canted spin 
arrangement at H < H c : 



Dk 2 + 2SA 1 



1 

2S 



(D(k - k c0 f + 2(H' C - H)) - S 2 {Ql z f sin 2 9 



(22) 



where k c0 and H' c are given by Eqs. ^ and (f5]), respectively, cos# = H/(2SA(1 — 1/25) — Suj a), 9 is the canting 
angle between spins and the field, and we have discarded terms proportional to A similar to the last one in Eq. (fTO|) 
which are of the next order in 1/S. The second bracket in Eq. (j22f vanishes at k = k c o when H = H' c . The last 
term in Eq. (|22|) gives negligibly small correction to H' c . The spectrum (|22|) is unstable at H > H' c . This result was 
obtained first in Ref.— . It was proposed there that a noncollinear phase arises at T = in the interval H' c < H < H c , 
where H c and H' c are given up to quantum 1/S corrections by Eqs. (|21[) and ([5]), respectively. Meantime we show 
below that the energy of the collinear phase at H = H' c is higher than that of the noncollinear one and the value 
of the critical field H' c is smaller (if it is nonzero at all) than that given by Eq. ([5|) and derived from the spin-wave 
analysis. 



IV. SPECTRUM RENORMALIZATION AT H > H c 



Corrections of the first order in 1/S to the spectrum stem from diagrams shown in Fig. [51 To calculate them one 
should take into account terms H.3 and H.4 in the Hamiltonian which have the form 



n 3 

7^4 




1 

8^ 



ki+k 2 +k 3 =0 

E 

ki+k 2 +k 3 +k 4 =0 



z a^_ 1 (a ] _ 2 + a 2 j a 3 , 
[4A + 2(Ji + J 3 - 2J 1+3 ) - (Qf + Qf + 4Q^ 3 )] al ia l 2 a z a A 



+ [QT-Qf -2A] atj 



a 2 a 3 + al 2 aL 3 



a 4 , 



(23) 



(24) 



where we drop index k in Eqs. (|23p and 

We have found in accordance with the conclusion of Refi£ that the loop diagram presented in Fig. [2]^b) is much 
smaller than the Hartree-Fock one shown in Fig. [21a). One obtains for the contribution to Q(u>, k) from the Hartree- 
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Fock diagram 



n(u, k) 



^q ~ £ q 
Se a 



e ^E 

q 

^Er ( E * E « + Bkfi q) ( J q - Jk +q + Q" - QtU) 



^E 



E q~ £ q A _ 3 2 



B k -V 



2Se„ ^ q 



Bq / gq-e q A _ 3 
e q V 25 +A 2^ 



(25a) 
(25b) 
(25c) 

(25d) 
(25e) 



Term (a) in Eq. (j2"5")l is equal in the first order in l/S to BkA/2 and it cancels the last term in Eq. (fT§]) . Term (|25b) 
does not renormalize the bare spectrum in the first order in l/S. Term (125b ) gives rise to the anisotropic corrections 
to the spectrum. 

It is convenient to use evident relations E^ = (E^ + Bk)/2 + (E^ — Bk)/2 and B^ = (E^ + B^)/2 — (E^ — B^)/2 
to extract combinations proportional to E^ — B^ and E^ + from terms (c)-(e) of Eq. (|25p . Terms proportional to 
-E-k — -Bk and E^ + B^ renormalize the first and the second brackets in Eq. (fT9)) . respectively. As we are interested 
in the phase transition to the noncollinear phase, we focus on renormalization of the second bracket in Eq. (|19p and 
consider only terms proportional to E^ + £?k- 

We assume for definiteness in the particular calculations below that J>j4>wo- As a result of simple but tedious 
calculations one has from Eq. f|25[) in the leading order of u>q/ J and A/ J for the important for us terms 



A 

—B k 

2 k 



(£k + Bk) (^k 2 loo [X cos(20 k + 2/3) cos 2/3 + Y cos 2 fa] + ^ E f 5 



X 



Y = 



D 



q ^q i 
— T y 1 - (cos q x - cos g z ) , 



165^0^ ^ Jo — Ja 
q M 



3 /2AS, A AST> 

2^v^ ln r(w 



(26) 
(27) 
(28) 



where /3 is the angle between magnetic field and an edge of the square (see Fig. [3]). The constant X should be 
calculated numerically because summation over large momenta is important in Eq. (|27[) and one cannot use Eqs. ([8]) 
and ([9]) for the dipolar tensor components. This calculation can be carried out using the dipolar sums computation 
technique (see, e.g., Ref^ and references therein) with the result X ss 0.0087 for exchange coupling between only 
nearest neighbor spins on the simple square lattice. 

Terms in Eq. ([26| proportional to X and Y stem from parts of the term (c) in Eq. (|25|) proportional to J q — Jk+q 
and <3 q z — <5 k + q , respectively. They make negligibly small correction to the constant D in the second bracket in 
Eq. (|19p . Meantime they are very important for our consideration because they fix the direction of the momentum 
at which the spectrum instability arises: k c o is directed so that this term to be minimal at fa = fa c0 ■ Notice that 
anisotropic terms are invariant according to replacement fa i— > fa ± 7r. Then k c0 and — k c0 are equivalent. One 
obtains after simple calculations that the anisotropic part in Eq. (|26|) has a minimum at fa o0 given by 



tan2</>k c0 = — 



sin 4/3 



r/A + 2cos 2 2/3' 



(29) 



where the solution should be taken with cos2</>k c0 < 0. The term Y/X plays in Eq. (|29|) only when j3 ~ ir/4 because 
X 3> Y. Simple calculation shows that the first and the second terms in denominator in Eq. (|29|) are equal at 
(3 = 7r/4 ± (5, where 



(30) 




« 1. 
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(a) < p < n/4 - 8 n/4 + 8 < p ^ 




FIG. 3: (Color online). It is assumed in this paper that the field H lies within the plane of the lattice and is directed by an 
angle to the square edge. This figure shows the direction of the vector k c o at which instability arises of the classical spin-wave 
spectrum at H < H c . Domains in the stripe phase are elongated perpendicularly to k c o. Slide (a) is for the case when /3 lays 
outside the narrow interval (7r/4 — S, n/4 + S) (cross-hatched area in the figure), where <J <§C 1 is given by Eq. (|30[1 . Slide (b) 
illustrates the rotation of k c o when the field rotates inside the interval (n/4 — S, n/4 + 5). The direction of the corresponding 
rotations are indicated by dashed arrows. We have k c o _L H when /3 = n/4. Notice that k c o and — k c o are equivalent. 



Then, when (3 lays outside the narrow interval (tt/4 — S, n/4 + S) the term with Y can be discarded in Eq. (f2U)) and 
one has for oV,, 



n/2-0, if/9e[0,7r/4-*), r . 

-P, if 13 e {n/4 + S,n/2}. 



Eq. (p{l"j) signifies that k c0 does not rotate with the field and is directed along the square edge as it is shown in Fig.[3{a). 
On the other hand, the term with Y comes into play when j3 lays in the interval (7r/4 — S, 7r/4 + S) and k c0 rotates in 
the same direction as the field H does so that k c0 turns out to be perpendicular to H when (3 — n/4 (see Fig. EJb)). 
The last term in Eq. (|26[) leads to quantum renormalization of the critical field H c which has the form 

nc = H?>-^(£z) (32) 



4nD \2AS 

when \n(D/AS) > 1, where H ( C Q) is given by Eq. (|21[) . It is seen from Eq. p2p that quantum fluctuations reduce the 
value of the saturation field. As a result we have for the spin-wave spectrum at k <C 1 

e 2 k » (Dk 2 + H) (D(k - k c0 ) 2 + A + 5(4> k )) , (33) 
A = H-Ho, (34) 



where H c is given by Eq. (|32| and <5(0k) stands for the small anisotropic terms discussed above. 
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V. THE NONCOLLINEAR PHASE 



We discuss in this section properties of the noncollinear phase. Spin arrangement is considered at H w H c and 
H < H c in subsections IV Al and IV Bl respectively. The spin- wave spectrum is discussed in subsection IV CI 



A. Transition from the collinear phase to the noncollinear one at H — H c 

It it shown in the previous section that the spectrum of the collinear phase becomes unstable when A < 0, 
i.e., when H < H c . The instability takes place at momentum k c o which absolute value is given by Eq. ^ and the 
direction is determined by small l/S'-corrections. This instability signifies a transition to a noncollinear phase. We 
discuss now the spin arrangement in this phase using the approach of Bose-Einstein condensation of magnons. 8 

There are two equivalent minima in the bare spectrum (|33|) at momenta ±k c o. The instability of the spectrum at 
k = ±k c o at H slightly smaller than H c means a "condensation" of magnons at states with these momenta. It is 
shown below that at H < H c the energy can be lower by varying the value of the momentum k c ||k c o at which the 
condensation takes place. Thus, k c depends on H so that k c = k C Q at H = H c and k c < k c o at H < H c . Then we 
imply below that k c ^ k C Q. 

At H < H c one should replace operator a± kc and a. j_ k as follows 



a±k c i-> a±k c + 



(35) 



where p ka is the density of condensed particles and (p± is a phase constant which should be found from the demand that 
linear terms disappear in the Hamiltonian arising after transformation (|35[) (see, e.g., RefJ^). Taking into account 
the equivalence of minima at +k c and — k c we assume in Eq. (|35[) that p kc = p- ka and similar equalities are implied 
below for harmonics of k c . One has from Eqs. (Ti"2"| and (|24p for the term in the Hamiltonian linear in 

V^PkT [e l ^E kc + e-^B kc + 3Ap kc (e^ - - l^J) a k c - ( 36 ) 

We obtain the term linear in ak c conjugating Eq. (|36|) . In order expression (|36p to be equal to zero at nonzero positive 
Pk c , <Px should be equal to ±7r/2. We choose plus sign below and one obtains that Eq. (f36|) is equal to zero when 



7T 

Pi 



Pk c 



2 ' 

E kc - B kc _ D(k c ~k c0 ) 2 +A ( 3? ) 
6A 6A 



It is seen that positive solution for p kc exists only if A < 0, i.e., at H < H c . Notice that the solution (|37[) corresponds 
also to a result of minimization with respect to (fx and pk c of the correction to the energy arising after transformation 
(|35|) and given by 

= 2 Pkc {E kc + B kc cos 2(pi) + 3Api c (1 - cos 2y»i). (38) 

Minimization of £\ with respect to k c gives k c — k C Q. The deviation of k c from k C Q arises after taking into account 
condensation of magnons at states characterized by harmonics of k c . 

It is seen from Eq. (|24|) that transformation (|35|) leads also to terms in the Hamiltonian linear in a±3k c and aj_ 3k . 
In order to cancel these terms one has to perform the following transformation of a±3k e and a± 3k similar to (|35[) 



a±3k c i-* a±3k c + e %V3 \/9Tp3k c 



' - ^al 3kc +e-^ 3 ^/mp~3 



l ±3k c + e V Jl P3k c - 

As a result the term in the Hamiltonian linear in arising after transformation (|39p has the form 

3/2 



(39) 



^ (x/pik7 (e ltp3 E 3kc + e-^B 3kc ) + 2iAp^ + Ap kc ^pikT (7e^ 3 - Se"*-)) a k- ( 40 ) 
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It is equal to zero if 

TT 

¥>3 = - 77, 

, 1 D(k c -k c0 ) 2 +A (41) 



3 D(3k c - k c0 ) 2 - 2D{k c - k c0 ) 2 - A ' 



Solution (|4ip corresponds also to the result of minimization with respect to tp 3 and p3k c of the correction to the energy 
arising after transformation (|39|) that has the form 

7^- = %P3kc( E 3k c + B 3kc cos2tp 3 ) + 8Ap^ c 2 y/p 3kc sintpi smip 3 + 2A(7 - 5 cos 2tp 3 )p kc p 3kc . (42) 

Condensation of magnons at states characterized by harmonics of k c makes it possible to lower the energy further 
by varying the absolute value of k c at any given H . To show this let us consider first terms in Eqs. (|38|) and ([42]) 
which are the only terms depending on k c in these expressions. Their explicit forms are 2p] SLa (D(k c — fc c o) 2 + A) and 
2p3k c (£ ) (3fc c — k c o) 2 + A), respectively. It seen that at any fixed pk c and p3k c a small deviation of k c from k c g (to 
be precise, small decreasing) raise E\ and lower £ 3 . The total energy lower in this case because £\ rises quadratically 
with (k c — k c o) and £ 3 lower linearly with [k c — k c o)- As a result minimization of the total energy correction E\ + £ 3 
with respect to k c gives in the leading order in A/Dk 2 

A 

c c0/ 



k ^^( 1 -M[DkT ] 1 (43) 



Notice also that k c should be parallel to k c o because the absolute value of negative anisotropic corrections to D 
discussed in the previous section have a maximum when k c ||k c o. It is seen from Eqs. (|41[) and (|43[) that p 3 k c <C Pk c if 
H « H c . We use this fact in Eqs. (j4H)) and (|4"2")) discarding terms of powers in p^ higher than 1/2 and 1, respectively, 
and omitting terms with higher harmonics in Eqs. (f31)|) , (|3"8|) , P0|) and P2"| . 

It is easy to show that the phase factor of (2n + l)-th harmonic of k c should be equal to tt/2 or —tt/2 in order 
the linear terms can vanish in the Hamiltonian. One has to take into account the fact that terms in the Hamiltonian 
linear in a(2n+i)k c and a[ 2 n+i)k stem from TL q with even q and that B(2 n +i)k c ~ -B(2n+i)k c ~ SA if n is not very 
large. Closer examination gives 

Van+i = (-l) n f • (44) 

For instance, it is seen from Eq. (|4"2"1) that and ip 3 should have different signs to minimize the second term. The 
validity of Eq. (f44|) for higher harmonics is not so evident for large enough |A|, at which densities of condensed 
particles at states characterized by neighboring higher harmonics differ several times only. Meantime particular 
numerical calculations show the validity of Eq. (|44|) in this case too. 

Terms in the Hamiltonian with odd number of operators a and which are proportional to Fourier components 
of the non-diagonal parts of the dipolar tensor Q xz (see, e.g., Eq. (f2"3")0 lead to condensation of magnons at states 
characterized by even harmonics of k c . It is easy to show that phase constants (p2 n should be equal to zero modulo 
tt in corresponding transformations similar to (|35[) and (|39|) . The density of condensed particles are very small at 
states characterized by even harmonics of k c . For instance, p2k c oc p\ (Q^k M) 2 ^ Pk e - Notice that P2nk c = when 
k c _L H because Q^nk = hi this case. The density of condensed particles in the state with zero momentum po is 
equal to zero at any k c direction because Qg z = 0. 

One can find the spin ordering in the noncollinear phase at H « H c using Eqs. (|37|) . (|4Tj) and (|43"]) . We obtain 
mean values of the transverse spins components (S£) w v / 2S'(ak + al k )/2 and (3%) ~ V2S(a^ — a^ k )/2i using the 
fact that (a n k c ) = e* v " y/^Prik c - Only even and odd harmonics contribute to (Sf) and (Sf), respectively, because 
<y32ri+i = ±7r/2 and (p2 n — mod tt and we have 



(Sf) = y/8SpV c cos(k c R0 + ^(-l)"J P(2 " +1)kc cos((2n + l)k c R0 . (45) 

\ n— 1 ' P^~° / 

Due to the smallness of densities of condensed particles at states with even harmonics we do not discuss (Sf ) here. As 
P(2n+i)k c ^ Pk c when H w H c , only the first few terms play in Eq. (|4"5)) at such H. We plot the value (Sf) / ^8Spk c 
in Fig.g]for H w iJ c and H = H c - l.5Dk 2 using Eqs. ([27]), (HTJ, (0S1), and (gSJ). It is seen that (5,f) has a sinusoidal 
profile at H ss iJ c as it is shown in Fig. [1] that is slightly distorted at smaller field due to the higher harmonics in 
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FIG. 4: (Color online). The value {Sf)/^8S^T C given by Eq. ((45} is shown in the noncollinear phase (i.e., at H < H c ) to 
demonstrate the field evolution of the domain pattern profile forming by off-plane spins components (see Fig. [TJ . It is seen that 
the profile is sinusoidal at H « H c that evolves upon the field decreasing into a meander-like and then into a step-like form. 
Notice that the period of the stripe structure rises upon the field decreasing as it is illustrated in Fig. [S] 



Eq. P5| . Higher harmonics in Eq. (f4"5)) become important at H < H c — 1.5Dk^ when they change considerably the 
profile of (Sf) (i.e., the domain pattern profile) which we discuss in the next subsection. 

Let us compare energies of the stripe phase and the collinear canted phase at H < H c . The last one can be found 
performing the transformation for operators ao and a} Q similar to (|35[) and (|39p and minimizing the corresponding 
expression for the energy correction. The result is 

£o_ (E -B )2 _ GPfc c 2 o-|AQ 2 U6) 
m AA AA { ' 

This expression gives the value of the energy correction higher than that of the stripe phase. For example, correspond- 
ing values for H = H' c are — 0.179A 2 /A and — 0.063A 2 /A for the stripe and the canted collinear phase, respectively, 
where H' c given by Eq. ([5]) is the field below which the classical spin-wave spectrum of the canted phase is stable. The 
difference between energies of the collinear canted phase and the stripe one decreases as the field reduces. 



B. H < H c 



We demonstrate in Appendix [S] that expression for the energy of the noncollinear phase is much more complicated at 
H < H c because one has to take into account corrections to the energy containing P(2«+i)k c with n > 1 for accurate 
finding of the stripe profile and k c . The corresponding values of P(2n+i)k c an d k c have to be found by numerical 
minimization of the energy. The number of harmonics of k c to be taken into account in particular calculations rises 
as the field decreases. For example, we find that one can use only five harmonics for H = H c — 2.5Dk^ whereas their 
number should be increased to 51 for H = H c — 16Dk^ . This circumstance leads to a restriction on the field interval 
which can be considered practically. The smallest field value we have managed to reach is H = H c — 16Dk^ Q . 

We plot (Sf )/\/8S'pk c in Fig. Q] for some H values using Eq. P5"|) and results of numerical calculation of P(2 n +i)k t .- 
We used Eqs. ([37]), fH]), (gSJ for H = H c -1.5Dk% and Eqs. jMJ) for H = R c -2.hDk 2 ^. It is seen that the sinusoidal 
profile evolves upon the field decreasing into a meander-like and then into a step-like form. Then, the domain wall 
density decreases from unity to a value much smaller than unity upon the field decreasing. 

Dependence of k c on H is illustrated in Fig. Notice that the period of the stripe structure given by 27r/fc c rises 
very fast upon the field decreasing. This finding is in qualitative agreement with results of Ref.— , where a very large 
period of the stripe pattern was found at H = when A 3> ujq. 
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Dkl, 

FIG. 5: Dependence of the absolute value of the momentum k c describing the stripe structure on the magnetic field found by 
numerical energy minimization taking into account 51 harmonics of k c . The inset shows the neighborhood of point H = H c , 
where the dependence of k c on H is described by approximate analytic expression (|43p (dashed line). 



C. Spin-wave spectrum in the noncollinear phase 



There are two types of corrections to the bilinear part of the Hamiltonian (|12|) arising after transformation (|35|) . 
The first one leads to the following simple rcnormalization of Ek and B^: 

E k ^E^ + 7Ap kc , S k h-> 5 k - 5 A Pkc . (47) 

This renormalization gives rise to an additional term equal to 2|A| in the second bracket in Eq. (|33|) and the spectrum 
takes the form 

e k « (Dk 2 +H) (D{k-k c ) 2 + \A\). (48) 

Higher harmonics contribute to renormalization of and B^ as well but their effect is small at |A| < Dk 2 because 
corresponding densities are much smaller than p^ c . 

Another sort of corrections to the bilinear part of the Hamiltonian include the so-called umklapp interaction having 
the form 



7 + •> / f f 

+ a k°-k±2k £ 



(49) 



There are also terms given by Eq. (l4"9")) with 2nk c and /o„k c instead of 2k c and pk c with n > 1 which can be neglected 
at |A| < Dk 2 . The umklapp interaction makes much more complicated the spin-wave spectrum analysis (see, e.g., 



RefJ£ for some detail) . Meantime it can be shown that corrections to the spectrum square from umklapp terms 
are negligible for momenta which absolute values are not very close to k c : D(k — k c ) 2 > |A|. Then the spectrum is 
given by Eq. (gSJ) in this case. Carrying out of much more complicated analysis for k m k c is out of the scope of the 
present work. 



VI. CONCLUSION 



In conclusion, we discuss 2D Heisenberg ferromagnet with dipolar forces and out-of-plane easy-axis one-ion 
anisotropy in in-plane magnetic field described by the Hamiltonian (P). We consider the classical spin- wave spec- 
trum assuming a collinear spin arrangement and find in accordance with previous results^ that the spectrum is 
unstable in a range of the field H' c < H < H c , where H c is the saturation field given by Eq. ([52^1 and H' c is given 
by Eq. ([5]). The instability arises at an incommensurate momentum k c o which direction is not fixed in the linear 
spin- wave approximation and which absolute value is given by Eq. ([5]). Thus, we conclude following Ref. 6 that the 
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collinear spin arrangement is unstable at H < H c with respect to a stripe domain formation with a large period (see 
Fig. P. 

To discuss properties of the non-collinear phase we consider first 1 / S corrections to the spectrum of the collinear 
phase at H > H c and find that they fix the direction of k c o so that it depends on the field direction as it is shown in 
Fig. [3l We discuss the noncollinear phase using the technique of the Bose-Einstein condensation of magnons suggested 
in Ref.— . Appearance of the stripe spin arrangement corresponds in this technique to a " condensation" of magnons 
at H < H c at states characterized by a momentum k c ||k c o and its harmonics (i.e., nk c , where n is integer). It is 
the unusual result of the present analysis that the vector k c varies with the field (remaining parallel to k c o) so that 
k c = k c o at H = H c and k c < k c o at H < H c . The origin of this behavior is in condensation of magnons at states 
characterized by harmonics of k c . 

We demonstrate that the energy of the stripe phase is lower than the energy of any collinear spin arrangement 
at H < H c but the difference between these energies decreases as the field reduces. We have not found any sign of 
transition to a collinear phase at H < H c . Then, H' c given by Eq. ([5]) is overestimated in the previous work£ in which 
it was obtained from the condition of the spin-wave spectrum stability in the collinear canted phase. It should be 
noted that this our finding is in qualitative agreement with results of Ref.~, where it was shown, in particular, that at 
H = the difference between energies of the stripe pattern and the collinear phase is very small in the limiting case 
of A ^> luq discussed in the present paper. 

Expression for values of spins components perpendicular to the film forming the domain pattern is given by Eq. (|45D . 
Coefficients /9„ kc and the value of k c are obtained by numerical minimization of the ground state energy. The 
transformation of the domain pattern upon the variation of the field is shown in Fig. 0J It is seen that the domain 
pattern has a sinusoid profile at H ss H c that transforms to a step-like profile upon the field decreasing so that the 
domain wall density decreases from unity to a value much smaller than unity. The dependence of k c on the field is 
shown in Fig. O It is seen that the period of the stripe pattern (27r/fc c ) rises very fast as the field decreases. This 
finding is also in qualitative agreement with results of Ref£, where it was shown that at H = the period of the 
stripe structure is very large at A ^> luq. 

We find that the spin-wave spectrum in the noncollinear phase is given by Eq. (j4"5|) for momenta D(k — k c ) 2 > |A| 
if |A| < Dk%. 
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APPENDIX A: ENERGY OF THE NONCOLLINEAR PHASE 

We present in this appendix some detail of the noncollinear phase energy calculation that was used in finding the 
value of k c drawn in Fig.[5]and the domain pattern profile presented in Fig.[4j It is convenient to express the correction 
to the energy arising after transformations ([33]) and (j3"9"|) (and similar ones for higher harmonics of k c ) in the following 
form: 

oo 

S = J2^n+l, (Al) 
n=0 

where £2n+i contains all terms depending on P(2n+i)k c an d not depending on P(2i+i)k c with i > n. In particular, one 
has for the first three terms in Eq. (|Aip 

E x = 2p kc (D(k c ~ k c0 ) 2 - |A|) + 6Apl c , (A2) 

E 3 = -8A P 3 k / c 2 ^r c +2AAp kc p 3kc +6Ap 2 3kc +2p 3kc (D(3k c ~k cQ ) 2 - |A|) , (A3) 
E5 = -24Ap kc% /p 3kc p 5kc + 24v4p 3kc ^/p kc/ 9 5kc + 24A(/9 kc + p 3 k c )p5k c + 6Apl kc 

+ 2p 5kc (D(5k c - k c0 ) 2 - |A|) , (A4) 

where we take into account Eq. (|4~4"1) for phase factors (fi2n+i- We can restrict ourself by some first terms in Eq. (|A1[) 
to perform the numerical minimization of E because the value of P(2n+i)k c decreases as n rises. Meantime, the number 
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of terms to be taken into account in Eq. (|A1[) rises as the field decreases. For example, we find that one can use only 
terms (|A2I) — (|A4[) for H = H c — 2.5Dk^ Q . The corresponding values of k c and densities are given by 

k c fa 0.832fc c0 , (A5) 

Pkc „ 1.13B, ./^Ir, 0.183, .[S « 0.0325, ./^» 0.0057, (A6) 

V P^c V V Pkc 

which, in particular, have been put into Eq. lf4"5]) to plot (Sf) /^/85pk c in Fig. [U We have obtained that one has to 
take into account 26 terms in Eq. (IA1|) for H = H c — 16Dk^ Q . 
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